Is Small Perfect? Size Limit to Defect Formation in Pyramidal Pt Nanocontacts 



V. RodrigucsiQ F. Sato^, D. S. Galvao\ and D. Ugarte^-^ 
^ Instituto de Fisica "Gleb Wataghin", Universidade Estadual de Campinas, 
Unicamp, CP. 6165, 13083-970, Campinas, Sao Paulo, Brazil and 
^Laboratorio Nacional de Luz Smcrotron, CP. 6192, 13084-971, Campinas , Sao Paulo, Brazil 

(Dated: February 5, 2008) 

We report high resolution transmission electron microscopy and ab initio calculation results for 
the defect formation in Pt nanocontacts (NCs). Our results show that there is a size limit to the 
existence of twins (extended structural defects). Defects are always present but blocked away from 
the tip axes. The twins may act as scattering plane, influencing contact electron transmission for 
Pt NC at room temperature and Ag/Au NC at low temperature. 
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The structural and mechanical properties of nanomet- 
ric wires represent a fundamental issue for the under- 
standing of different phenomena such as friction, frac- 
ture, adhesion, etc. [1]. There is a renewed interest 
in such systems motivated in part by the growing de- 
mand on miniaturization, new functionalities and less 
power consuming systems. Components such as diodes 
[2], switches and electronic mixers [4 have been built 
with simple molecules as active units [5\ . In spite of these 
technological advances, important issues on device inte- 
gration, how to connect them in a stable and reproducible 
way, are still open [6 , 7J . In order to build functional nan- 
odevices, electrical contacts of nanometric dimensions are 
needed. The contact atomic arrangement will probably 
determine the electronic structure and, in consequence, 
the coupling between leads and the active part (for exam- 
ple, molecules) of the device [7J[S1[3- this context, it is 
necessary to understand the influence of contact proper- 
ties in the device characteristics and also the role played 
by structural defects since they can compromise the de- 
vice functionality and reliability. Unfortunately, we can 
not control yet the nanocontact (NC) fabrication and, its 
characterization is rarely possible; consequently, control 
and reproducibility are still very difficult, even for the 
simplest cases, such as the hydrogen molecule jTOj . 

Metal conductor NC studies have been usually done 
in experiments where the NCs are generated by mechan- 
ical elongation. In order to get a deeper understand- 
ing on the NC properties, its mechanical and atomic 
structure should be better comprehended. For example, 
for atomic size NC generated by mechanical stretching, 
we can expect that the metal deformation must exhibit 
stages where structural defects are formed to acomodate 
the elongation strain. It has been already shown that 
NCs obtained by mechanical elongation exhibit distinct 
elastic strain stages, followed by sharply deflned yield- 
ing, originated from structural reorganization [TI ITTIIT^ . 
In this way, we can not properly address the deformation 
of nanosystems using continuum-based theories, a micro- 
scopic analysis taking into account the atomic structure 
is necessary. Today, most of atomistic process informa- 



tion comes from molecular dynamics simulations where 
the elongation velocities (m/s) are orders of magnitude 
larger than in experiments (nm/s) [U [13l [TH [15]. The 
plastic deformation is associated with collective slips of 
entire atomic planes or order-disorder transitions |12l 113] . 
It has been suggested that at very small scales (10 nm) 
dislocations should be completely suppressed because the 
involved stress are comparable to the intrinsic lattice 
strength. In analogy, a self puriflcation effect has been 
recently reported for semiconductor particles: defect for- 
mation energy increases as the nanocrystal size decreases, 
rendering difficult the doping of smaller systems [ini [TT] . 

In this sense some fundamental questions must be ad- 
dressed. Is small perfect? Can we generate defect free 
metallic NCs? Is there a size limit to the existence of 
such structures? In this letter we present results from dy- 
namic high resolution transmission electron microscopy 
(HRTEM) experiments of the atomic arrangement of Pt 
NCs under mechanical deformation at room tempera- 
ture. We have observed that there is a size limit where 
the twins (extended structural defects) can exist, below 
this limit the defects are always present but blocked few 
atomic planes away from the atomically sharp Pt NC. 
In order to quantitatively estimate the size limit where 
defect free contacts can exist, we have estimated the ener- 
getics associated with the twin generation through ab ini- 
tio calculations using the well-known and tested SIESTA 
code [TH]. 

The atomic size metal NCs were generated in situ in a 
HRTEM (JEM 3010 URP, 300 kV, 0.17 nm point resoL, 
operating at LME/LNLS, Campinas, Brazil). In this 
methodology, the microscope electron beam is focused 
to generate holes in a self-supported thin metal film un- 
til neck is formed between them [T^] . The NCs (1-2 nm in 
diameter) spontaneously become thinner due to displace- 
ments of the apexes, reach the size of a few atoms and 
finally break. During this process, we have no control on 
the direction or speed of the elongation. Real time evo- 
lutions were recorded using a high sensitivity TV camera 
(Gatan 622SC, 30 frame/s). 

For the theoretical analysis, we have used ab initio den- 
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sity functional theory (DFT) in the framework of the 
local density approximation (LDA) [18| . In order to war- 
ranty high precision results, we have used double- basis 
set plus polarization functions (relativistic calculations) 
and norm-conserving pseudopotentials built on TrouUier 
Martins scheme [5D]. We first tested the parameters used 
on SIESTA code for Pt perfect crystalline bulk struc- 
tures, relaxing both cell lattice vectors and atomic posi- 
tions. The obtained values are in good agreement with 
the experimental ones (3.924 and 3.875 A , for the ex- 
perimental |21j and SIESTA values, respectively). Once 
established the pseudopotential reliability for Pt struc- 
tures, we proceeded to analyze the energetics and relative 
stability of the structures with defects (twins) . 

The analysis of numerous HRTEM images always dis- 
play defect free crystalline Pt nanoconstrictions com- 
posed of few atomic planes. When present, defects arc 
observed away from the tip. As an illustrative exam- 
ple, figure [l] shows a nanostructure where the lower apex 
is well aligned with the HRTEM beam, thus providing 
atomic and time resolved images of a pyramidal tip con- 
tact. Its axis is along the [111] crystalline direction (here- 
after called [111] tip) and that ends with just only one 
atom. In the lower apex, the horizontal lattice periodic- 
ity is 0.24 ± 0.01 nm (Pt (111) spacing is 0.227 nm) and 
the angular relations confirm that the pyramid is defined 
by a (111) and (100) facet at left and right, respectively. 
The apex vertical movement generates the NC thinning, 
but a small shear movement (upper apex is gliding to 
the right) has induced the formation of a kink on the left 
apex side in (a) . Further deformation generates an atom- 
ically sharp tip formed by 7 atomic layers (c); the arrow 
indicates the formation of a twin defect at the 5th atomic 
layer (counting from the tip); after rupture (f), the twin 
defect is annihilated and the tip reorganizes itself to form 
a truncated pyramid. See supplementary material in |22j . 




FIG. 1: Sequence of images showing the elongation and rup- 
ture of a Pt NC under tensile stress with a minor shear com- 
ponent (observation axis [110]). 



From figure [T] we can also observe the presence of pla- 



nar defects, more precisely a twin defect. This can be 
more precisely identified by considering that a twin de- 
fect is a fault in the ABC planes sequence forming a 
mirror: ABC B AC. On the experimental image it can be 
recognized by the angle change in the pyramidal facets. 
Above the twin, the left facet of the tip is a (111) plane 
and a (100) one on the right; bellow the twin it is a (111) 
plane on the right and (100) on the left (see figure |2]). 

One very important result is that during the NC for- 
mation no twin was observed very close to the tip. It has 
only been observed at the 4th layer (figure [ij^d)) and at 
the 5th layer (figure [ije)). These results raise the im- 
portant question whether this represents the transition 
size between the macroscopic plasticity mechanism and 
the nanometric scale where extended defects can not be 
sustained anywhere [T]. In order to better address this 
question, we have analyzed the energetics associated with 
this defect formation. We have calculated the structure 
relaxation and the total energy of the structures shown in 
figure [3j a)- (e). As the experimentally observed tip (fig- 
ure [ij has seven layers, we have investigated all the five 
possible positions for the twin defect, as shown in figure 

EI 




FIG. 2: (a) NC with a twin defect formed in its 5th layer 
counting from the top. (b) Schematic draw of the tip in (a) 
showing the usual ABC notation for the face-centered-cubic 
packing. 



In order to model the tip based on our experimental 
data we should consider that the HRTEM images pro- 
vide a bi-dimensional projection of the atomic arrange- 
ments, but tree-dimensional data is needed to build in- 
put models for the calculations. The apex morphology 
is determined by the crystal faceting |TT], which can be 
easily obtained using the Wulff construction rules [53] ■ 
For the Pt structures in our case (figure [ij only (111) 
and (100) facets needed to be considered. However, ad- 
ditional factors must be taken into account in order to 
build an atomically sharp [111] tip. For example, the 
atom at the extreme of the tip is located at the center of 
a regular triangle (2nd layer), which by itself is located 
over an hexagonal structure (3rd layer, figure [sjj a)). Al- 
though in this configuration the atom that finishes the 
tip is not in the expected crystallographic position, this 
is the unique way to construct the sharp pyramid. The 
triangular shape of the 2nd layer fixes the apex [100] facet 
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width to be two-atom-row in size (figure|3ja), named AT2 
apex). To generate the twinned apexes, we must consider 
two crystaUine domains; a base and a tip. A twin can be 
also seen as a 60° rotation of one of the domains along 
the tip axis; in this way, the (111) (or (100)) facets of the 
tip become (100) (or (111)) facets at the base (see figures 
[2] and [3]). Although the domains are rotated, they must 
have exactly the same cross-section at the twin. Briefly, 
the faceting pattern is determined by two geometrical 
constrains: a) an atomically sharp apex imposing that 
the tip domain must be identical to the equivalent re- 
gion of the AT2 apex; b) the twin position, because the 
tip (100) (or (111)) facet width at the twin position will 
determine the width of the (111) (or (100)) facet at the 
base. As a consequence, these apexes will not always 
follow rigorously the Wulff surface balance and also, the 
total number of atoms in the tips will not be constant 
(see Tabic 
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FIG. 3: Schematic draw of possible atomic arrangement of 
twinned [111] Pt NCs. The schemes in (a)-(e) are organized 
as a function of the twin height position at the apex, 2nd, 
3rd, 4th, 5th and 6th layer respectively. Lines indicate the 
atomic layer and the facet borders. 



Based on the procedure discussed above, we generated 
the 3D structures for the calculations. The initial Pt 
atomic positions followed a geometrically (or twinned) 
Pt lattice as schematically shown in figure [3] and then 
we have carried out fully geometrical optimizations. In 
Table|l]we present the relative energy formation: Egtruc— 
NEa, where E^truc is the energy of the configuration, N 
is the number of the atoms in the structure, and Ea the 
energy of one isolated atom. The results are relative to 
the lowest energy configuration (AT4) . 

Theses results can be explained in terms of faceting and 
stress release energy generated by the twin defect [T1E5]. 
AT4 is the favored configuration because it allows, at the 
same time, one-atom tip and an optimized base surface 
(figure [3|c)). Moving the defect one or two layers up 
(AT3 and AT2, respectively) cost a large amount of en- 
ergy (47.02 and 46.57 eV, respectively). These structures 
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Au 
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of atoms 
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(eV) 


2(AT2) 










3(AT3) 


99 


47.02 


32.20 


34.59 


4(AT4) 


105 


0.00 


0.00 


0.00 


5(AT5) 


105 


1.35 


0.51 


0.36 


6(AT6) 


102 


26.53 


17.90 


18.99 



TABLE I: DFT formation energy estimated for twinned Pt, 
Ag, and Au NCs, with relation to the AT4 tip. See figure 3 



present imperfect faceting but the structural instability 
comes mainly from the loss of 6 atoms. In contrast, mov- 
ing the defect one layer down (ATS) changes the struc- 
tural energy by only 1.35 eV. However, moving one layer 
further (AT6) requires a large amount of energy (26.53 
eV), but considerably less than in the AT2 and AT3 cases. 
Also, the transition from AT5 to AT6 involves only the 
loss of 3 atoms. 

The AT5 apex represents the first available structure 
to absorb accumulated stress; this is in excellent agree- 
ment with our time resolved HRTEM observations which 
show that AT5 (figure [T]^d)) is formed much before the 
NC thinning and rupture. As for the possible structural 
transition from AT4 to AT5 apexes, it is important to 
note that we are dealing with very small systems with a 
large ratio surface/ volume atoms. The twin can be con- 
sidered as a rotation between two configurations, this is 
certainly a very low cost process and probably with a 
low energy barrier for the size of the structures consid- 
ered here, as described in the quasi-melting of small Au 
clusters [H] . This interpretation can be better evidenced 
from video in supplementary material in [221 . 

In order to verify whether the obtained results for Pt 
could be extrapolated to other metals, we repeated the 
calculations for Au and Ag, those of most used metals in 
nanocontacts. The results are displayed in Table [l] and, 
the conclusions are identical, suggesting that these are 
general features of fee metallic nanocontacts. However, 
it must be remarked that generated Pt twins are blocked 
by a quite high energy barrier [? ? ], much higher than 
themal energy at room temperature. So, calculations ne- 
gleting temperature effects, such as the one presented 
here, can not be considered representative and a reliable 
discription of room temperature Pt NC experiments. On 
the other side, energy barrier are much lower for Ag and 
Au [? ? ]; this fact can account for the HRTEM results 
which show that defects are quickly annhilated at room 
temperature [HI [25]. But, it must be expected that the 
predicted Ag or Au size limit for twin defect formation 
should be observed in low temperature experiments 

In summary, we have shown that the elongation of Pt 
NCs induces the formation of twins located a few planes 
away from the apex, where the atomic arrangement of the 
atomically sharp Pt pyramid remains close to the ideal 
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one. This means that it is possible to produce one atom 
nanocontacts with a well defined structure. Many models 
in molecular electronics have assumed perfect tips, which 
has been object of criticism as being unrealistic. Our re- 
sults showed that this is in fact a good approximation. 
On the other hand, we always will have the presence of 
extended defects at the apex that would work as a scat- 
tering plane, thus it is expected that full electron trans- 
mission coefficients are not longer possible, at least for Pt 
NCs. These aspects should be taken into account when 
modeling the typical two point electrical measurement 
on molecular electronic devices. Similar results were also 
obtained for Au and Ag, suggesting that this is a gen- 
eral behavior of face-centered-cubic metals. Nevertheless 
they should be only observed at lower temperates (few 
Kelvin). It seems that 4 and 5 atomic planes is the limit 
size where twin defects can exist in Pt nanotips, thus rep- 
resenting the frontier to macroscopic plasticity regime for 
this metal at room temperatures. 

This work was supported in part by IMMP/MCT, 
IN/MCT, THEO-NANO, LNLS, CNPq and FAPESP. 
The authors acknowledge P.C. Silva and J. Bettini for 
assistance during the sample preparation and data treat- 
ment. 



[15] A. Nakamura, M. Brandbyge, L. B. Hansen, and K. W. 

Jacobsen, Phys. Rev. Lett. 82, 1538 (1999). 
[16] G. M. Dalpian and J. R. Chelikowsky, Phys. Rev. Lett. 

96, 226802 (2006). 
[17] S. C. Erwin, L. Zu, M. I. Haftel, A. L. Efros, T. A. 

Kennedy, and D. J. Norris, Nature 436, 91 (2005). 
[18] J. M. Soler, E. Artacho, J. D. Gale, A. Garcia, 

J. Junquera, P. Ordejon, and D. Sanchez- Portal, J. 

Phys.: Cond. Matter 14, 2745 (2002), see also 

http://www.uam.es/siesta. 
[19] Y. Kondo and K. Takayanagi, Phys. Rev. Lett. 79, 3455 

(1997). 

[20] N. Troullier and J. L. Martins, Phys. Rev. B 43, 1993 
(1991). 

[21] F. Cleri and V. Rosato, Phys. Rev. B 48, 22 (1993). 

[22] See EPAPS document No. X-XXXXX-XX-XXXXXX for 
the movie of the HRTEM experiment showing the evolu- 
tion of a Pt nanowire. A direct link to this document may 
be found in the online articles HTML reference section. 
The document may also be reached via the EPAPS home- 
page (http://www.aip.org/pubservs/epaps.html) or from 
ftp.aip.org in the directory /epaps/. See the EPAPS 
homepage for more information. 

[23] L. D. Marks, Rep. Prog. Phys. 57, 603 (1994). 

[24] W. Krakow, M. Jose-Yacaman, and J. L. Aragon, Phys. 
Rev. B 49, 10591 (1994). 

[25] V. Rodrigues, J. Bettini, A. R. Rocha, L. G. C. Rego, 
and D. Ugarte, Phys. Rev. B 65, 153402 (2002). 



* Electronic address: 'varlei@ifi.unic amp.br| 
[1] U. Landman, W. D. Luedtke, N. A. Burnham, and R. J. 

Colton, Science 248, 454 (1990). 
[2] R. M. Metzger and M. P. Cava, Molecular Electronics: 

Science and Technology (Annals of New York Academy 

of Sciences, 1998), pp. 95-115. 
[3] C. P. Collier, E. W. Wong, M. Belohradskyacute;, F. M. 

Raymo, J. F. Stoddart, P. J. Kuekes, R. S. Williams, and 

J. R. Heath, Science 285, 391 (1999). 
[4] J. Chen, M. A. Reed, A. M. Rawlett, and J. M. Tour, 

Science 286, 1550 (1999). 
[5] M. A. Reed, C. Zhou, C. J. MuUer, T. P. Burgin, and 

J. M. Tour, Science 278, 252 (1997). 
[6] H. B. Akkermanv, P. W. M. Blom, D. M. de Leeuw, and 

B. de Boer, Nature 441, 69 (2006). 
[7] J. Repp, G. Meyer, S. Paavilainen, F. E. Olsson, and 

M. Persson, Science 312, 1196 (2006). 
[8] S.-H. Ke, H. U. Baranger, and W. Yang, J. Chem. Phys. 

123, 114701 (2005). 
[9] N. Agrait, A. L. Yeyati, and J. M. van Ruitenbeek, 

Physics Reports 377, 81 (2003). 
[10] R. H. M. Smit, Y. Noat, C. Untiedt, N. D. Lang, M. C. 

van Hemert, and J. M. van Ruitenbeek, Nature 419, 906 

(2002). 

[11] V. Rodrigues, T. Fuhrer, and D. Ugarte, Phys. Rev. Lett. 

85, 4124 (2000). 
[12] G. Rubio, N. Agrait, and S. Vieira, Phys. Rev. Lett. 76, 

2302 (1996). 

[13] G. Rubio-BoUinger, S. R. Bahn, N. Agrait, K. W. Jacob- 
sen, and S. Vieha, Phys. Rev. Lett. 87, 026101 (2001). 

[14] A. M. Bratkovsky, A. P. Sutton, and T. N. Todorov, 
Phys. Rev. B 52, 5036 (1995). 



